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ABSTRACT 

We report on observations of a gamma-ray burst (GRB 061126) with an extremely bright (R « 12 
mag at peak) early-time optical afterglow. The optical afterglow is already fading as a power law 
22 seconds after the trigger, with no detectable prompt contribution in our first exposure, which was 
coincident with a large prompt-emission gamma-ray pulse. The optical-infrared photometric spectral 
energy distribution is an excellent fit to a power law, but it exhibits a moderate red-to-blue evolution 
in the spectral index at about 500 s after the burst. This color change is contemporaneous with a 
switch from a relatively fast decay to slower decay. The rapidly decaying early afterglow is broadly 
consistent with synchrotron emission from a reverse shock, but a bright forward-shock component 
predicted by the intermediate- to late-time X-ray observations under the assumptions of standard 
afterglow models is not observed. Indeed, despite its remarkable early-time brightness, this burst 
would qualify as a dark burst at later times on the basis of its nearly flat optical-to- X-ray spectral 
index. Our photometric spectral energy distribution provides no evidence of host-galaxy extinction, 
requiring either large quantities of grey dust in the host system (at redshift 1.1588±0.0006, based upon 
our late-time Keck spectroscopy) or separate physical origins for the X-ray and optical afterglows. 
Subject headings: gamma rays: bursts — gamma-ray bursts: individual: 061126 



1. INTRODUCTION 

While the study of the early-time X-ray afterglows of 
gamma-ray bursts (GRBs) has seen eno rmous strides 
since the launch of the Swift satellite (jGehrels et alJ 
2004), progress in the understanding of the longer- 
wavelength emission has been somewhat more measured. 
This is unfortunate, as a complete understanding of 
Swift afterglows can only come from a combined broad- 
band picture that allows us to systematically investi- 
gate whether the peculiarities seen in X-ray data carry 
over into the optical domain. Many of the same ques- 
tions raised by recent X-ray results can also be asked 
about the optical. Is there a prompt optical compo- 
nent, anal ogous to the steeply de caying component seen 
in X-rays (Barthelmy et al. 2005)? Does the optical light 
curve show unusual features suggestive of energy injec- 
tion, s uch as the nearly u biquitous X-ray shallow-decay 
phase (jNousek et al.l2006h ? Are there achromatic optical 
and X-ray breaks? Do the optical and X-ray afterglows 
even have a common origin at early times? 

Previous studies have provided important hints. Most 
observations have been interpreted to support the con- 
sensus picture of synchrotron emission originating from 
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a forward shock as it sweeps through the inter stel- 
lar medium (e.g., iDai fc Lulfl99l IVrba et alJ l2000h. or 
less commonly through a stell ar wind (|Price et al.l 1 20021 : 
iNvsewander et aT1l2006l ) — see lChevalierl (12007) for a re- 



view. In a smaller number of cases (lAkerlof et aUllMit 
Li et all l2003at iKobavashi fcZh ang 200l IShao k Pail 



20051 ). very early-time data have provided tentative evi 



dence for an additional emission component originating 
from the reverse shock as it travels backward into the 
shocked material in the frame of the forward shock. Re- 
cently, studies of early-time light curves have also shown 
evidence of significant delay betwe en the onset of th e 
prompt emission and the afterglow (Rvkoff e t al.ll200"4T ). 
and in at least one case complicated energy injection ac- 
tivity as late as nearly an hour after the gamma-ray burst 
(|Wozniak et al.ll2006l ). long after the gamma-ray emis- 
sion has faded away. However, simultaneous, correlated 
optical and gamma-ray emission has also been reported 
(iBlake et alj|2005t IVestrand et al JI20051 12006t IVost et all 
120071) for some events. In the Swift era, comparison to 
the very ea rly X-ray afterglow h as also been of great in- 
terest (e.g., lQuimbv et al.ll2006D . 

Most interpretations of the early afterglow have been 
based on unfiltered observations, or observations in a sin- 
gle filter. Without information about the frequency do- 
main, the reported early-time behaviors discussed above 
are difficult to definitively associate with any single phys- 
ical interpretation. Fortunately, the increasing number of 
fast-responding robotic ground-based observatories, the 
maturation of existing ones, and the rapid-response ca- 
pabilities of Swift are beginning to address this observa- 
tional gap. 

In the following discussion, we report on one of the 
brightest bursts of the Swift era, GRB 061126. The 
breadth and rapidity of the ground-based response to 
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this burst were remarkable, including unfiltered detection 
during the prompt emission and multi-color simultaneous 
detections in filters from U through K s (ranging a full 
decade in frequency) starting less than one minute after 
the burst trigger. This data set provides the opportunity 
to examine in unprecedented detail the time-dependent 
color properties of an early GRB afterglow. 

In Sj2] we present our observations from infrared (IR) 
through gamma rays of the early afterglow and our late- 
time Keck spectrum of the host, establishing the prob- 
able redshift of this system to be z = 1.1588. In 33.11 
we examine the properties of the prompt emission, and 
show that the high-energy and optical emission are obser- 
vationally uncorrelated temporally or spectrally, even at 
very early times. In S j3.2f43.3l we examine the properties 
of the optical-IR light curve, and provide evidence for a 
red-to-blue change in the spectral index of Af3 rj 0.3 at 
early times. We investigate the X-ray behavior in 34.11 
and show that no standard adiabatic model can fully ex- 
plain the behavior seen by the X-Ray Telescope (XRT). 
Finally, while in i j4.2H4.3l we show that the earliest after- 
glow appears reasonably fit by a reverse shock and the 
later afterglow by a forward shock based on the optical 
data alone, in 34 .41 we demonstrate that an extrapolation 
of the X-ray spectrum overpredicts the contemporaneous 
optical flux by a factor of 5-20. We demonstrate using 
the optical-IR spectral energy distribution (SED) that 
this discrepancy cannot be due to any known dust ex- 
tinction law. Unless we appeal to large quantities of grey 
dust, a possibility we discuss in §5.1) we argue in §5.2115.31 
that the X-ray and optical afterglow emission from this 
burst have separate physical origins. 

2. OBSERVATIONS 
2.1. Swift BAT and XRT 

At 08:47:56 on 26 November 2006 (UT dates are used 
throughout this paper), the Swift Burst Alert Telescope 
(BAT) triggered and located GRB 061126. Unfortu- 
nately, due to an Earth- limb constraint, Swift was unable 
to slew promptly to the target for 23 minutes and could 
not begin observations with the XRT or the Ultravio- 
let/Optical Telescope (UVOT) before that time. After 
23 minutes, Swift slewed to th e burst position and de - 
tected a fading X-ray afterglow (Sba rufatti et al.ir2006al) . 

We download the Swift BAT and XRT data from the 
Swift Archive 8 . The XRT data are processed with ver- 
sion 0.10.3 of the xrtpipeline reduction script from the 
HEAsoft 6.0.6 9 software release. We employ the latest 
(19 December 2006) XRT and BAT calibration files. We 
establish the BAT energy scale and mask weighting by 
running the bateconvert and batmaskwtevt tasks, also 
from the HEAsoft 6.0.6 software release. BAT spec- 
tra and light curves are extracted with the batbinevt 
task, and response matrices are produced by running 
batdrmgen. We apply the systematic error corrections 
to the low-energy BAT spectral data as advised by the 
BAT Digest website 10 . The spectral normalizations are 
corrected for satellite slews using the batupdatephakw 
task. 

8 ftp://legacy.gsfc.nasa.gov/swift/data . 

9 http://heasarc.gsfc.nasa.gov/docs/software/lheasoft/ . 

10 http: //swift . gsfc . nasa.gov /docs/ swift / analysis/bat_digest . html 



TABLE 1 

RAPTOR Observations of GRB 061126 





Filter 


Exp. Time 


Magnitude 13 


Flux b 


GO 




(s) 






23.4 


clear 


5.0 


12.26 ± 0.01 c 


38725.8 ±355.0 


32.3 


clear 


5.0 


12.66 ±0.02 


26791.7 ±489.0 


41.1 


clear 


5.0 


13.08 ±0.03 


18197.0 ±495.9 



Note. — Only the first three observations are given here. 
Additional data are available in the online supplement, or 
by request. 

a Exposure mid-time, measured from the Swift trigger 
(UT 08:47:56. 4). b Observed value; not corrected for Galac- 
tic extinction. Contemporaneous with strong GRB pulse; 
point not used in modeling. 

The reduction of XRT data from cleaned event lists 
output by xrtpipeline to science- ready light curves 
and sp ectra is described in detail in iButler fc Kocevs"kl 
(2007). The XRT, BAT, and RHESSI (Reuven Ramaty 
High Energy Solar Spectroscopic Imager) data are fit us- 
ing ISIS 11 . The XRT light curve is converted to normal- 
ized 100 keV flux using a scaling of 7.5 /zJy/cps. 

2.2. RHESSI 

RHESSI (|Lin et alj I2002D is a dedicated solar ob- 
servatory which uses nine germanium detectors to im- 
age the Sun at hard X-ray to gamma-ray energies (3 
keV - 17 MeV). These detectors are unshielded and 
therefore frequently detect emission from off-axis GRBs. 
GRB 061126 was detected by RHESSI, which with its 
large spectral range allows us to complete the high- 
energy spectrum of this event. 

To model the RHESSI response to off-axis photons, we 
have used Monte Carlo simulations and a detailed mass 
model. Since RHESSI rotates about its axis with a 4 s pe- 
riod, we have generated azimuthally averaged responses 
spaced 15° apart in polar angle. These responses are 
two-dimensional matrices of effective area: input photon 
energy vs. detected count energy bins. At present both 
energy axes are 64 logarithmic bins from 10 keV to 10 
MeV. We generate the response matrices with MGEANT 
(jSturner et al.! 2000) simulations: each response requires 
64 simulations of a monoenergetic input spectrum, one 
for each photon energy bin. For an individual GRB, we 
generate and subtract a background count spectrum us- 
ing data intervals before and after the burst. We generate 
a burst-specific response matrix with a weighted average 
of the two adjacent 15° responses. Convolving a spectral 
model with the response yields a model count spectrum 
for fitting to the GRB data. 

2.3. RAPTOR 

The RAPTOR ( Rapid Telescopes for Optical Re- 
sponse) experiment (|Vestrand et al.ll2002D consists of a 
series of small telescopes used to conduct transient sur- 
veys as well as perform rapid followup of GRBs and other 
events. One of these telescopes, RAPTOR-S, is a 0.4 m 
fully autonomous robotic telescope, typically operated at 
focal ratio //5. It is equipped with a lk x lk pixel CCD 
camera employing a back-illuminated Marconi CCD47- 
10 chip with 13 /xm pixels. The telescope is owned by Los 

11 http://space.mit.edu/CXC/ISIS/ . 
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TABLE 2 

PAIRITEL Observations of GRB 061126 



At a 


Filter 


Exp. Time 


Magnitude 13 


Flux b 


00 




00 




64.8 


J 


4.3 


11.72 ±0.02 


32658.7 ± 507.4 


262.0 


H 


4.0 


12.89 ±0.05 


7177.9 ±304.1 


136.7 


K S 


4.5 


11.01 ±0.05 


26278.5 ± 1298.0 



Note. — Only the first usable exposure in each band 
is given. The full table of photometry is available upon 
request, or in our online supplement. 

a Exposure mid-time, measured from the Swift trigger 
(UT 08:47:56.4) . b Observed value; not corrected for Galac- 
tic extinction. 

Alamos National Laboratory and located at the Fenton 
Hill Observatory (106.67° W, 35.88° N) at an altitude of 
^2500 m in the Jemez Mountains of New Mexico. 

RAPTOR-S responded to the Swift trigger at 
08:48:17.29, or 20.87 s after the trigger and 4.3 s af- 
ter receiving the GCN (Gamma-ray Burst Coordination 
Network) packet. The telescope took a series of nine un- 
filtered 5 s exposures (the first two of which occurred 
while detectable gamma-ray emission was still ongoing) , 
followed by a series of 10 s and 30 s exposures. The 
optical transient is detected in all these frames. Pre- 
liminary photometric calibration was performed using 
the i?-band magnitudes from the USNO (United States 
Naval Observatory) B1.0 catalog. However, for consis- 
tency with the unfiltered KAIT (Katzman Automatic 
Imaging Telescope) observations, which were calibrated 
using the more precise SDSS (Sloan Digital Sky Survey) 
measurements, we subtract a constant offset of 0.16 mag 
post-calibration. The RAPTOR photometry (not includ- 
ing this final offset) is given in Table [TJ 

2.4. PAIRITEL 

Starting in 2003, we began to automate the 1.3-m Pe- 
ters Telescope, fo rmerly used for the T wo Micron All Sky 
Survey (2MASS; ISkrutskie et ai]|2006l ). on Mt. Hopkins, 
Arizona. The telescope was re-outfitted with the south- 
ern 2MASS camera, and all of t he control and data acqui- 
sition systems were rewritten (|Bloom et al.l 12006). The 
Peters Automated Infrared Imaging Telescope (PAIRI- 
TEL) has been obtaining useful IR observations of GRBs 
since 2004. 

At 08:48:18 (f = 22 s), PAIRITEL was triggered with 
the GCN BAT position notice of GRB 061126. The 
autonomous slew of the telescope and dome began at 
08:48:22 and ended at 08:48:47; the slew time was short 
since we had been observing M82 (23.9° to the east of 
the GRB) immediately prior to the GRB. After an initial 
reset of the camera, the first 7.8 s images in J, H, and K s 
bands were obtained starting at 08:48:54.35 (t = 58 s). 
We continued with a dense sampling of observations over 
the next three hours as well as several hours of imaging 
the following night. 

Reductions of the individual images were performed 
using a set of customized scripts written in pyraf and 
Python. The afterglow was well detected with signal-to- 
noise ratio (S/N) > 10 in individual images for the first 
20 min of observations (Figure [1]). In fact, the H- and 
K s -band fluxes of the afterglow are so bright in the first 
few minutes that the pixel responses were in the non- 



TABLE 3 

NMSU 1 m Observations of GRB 061126 



At a 


Filter Exp. Time 


Magnitude 13 


Flux b 


(b) 




00 




(/■Jy) 


52.3 


I 


10.0 


12.53 ±0.00 


23746.5 ± 65.5 


96.8 


R 


10.0 


14.22 ±0.01 


6356.2 ± 35.0 


148.8 


V 


20.0 


15.32 ±0.01 


2737.8 ± 17.6 


213.2 


B 


40.0 


16.43 ± 0.01 


1108.7 ±9.2 


302.9 


U 


60.0 


16.67 ±0.03 


411.9 ±9.4 


Note 




Only the first filter cycle is 


given. Addi- 



tional data are available in the online supplement, or by 
request. 

a Exposure mid-time, measured from the Swift trig- 
ger (UT 08:47:56. 4). b Observed value; not corrected for 
Galactic extinction. 

linear regime. Unfortunately, the cloud cover in Arizona 
was highly variable during the first 30 min of GRB 061 126 
observations, leading to variable transmission on 10 s 
timescales. 12 As such, in our analysis we refit the zero- 
point in every individual exposure to the 2MASS cata- 
log. The typical root-mean-square (rms) uncertainties in 
the zeropoints are 2-3%. Given the large variations in 
the sub-pixel response function for NICMOS3 arrays, we 
have found that aperture photometry on individual ex- 
posures suffers a roughly 3% systematic uncertainty from 
image to image (Blake et al. 2007, in preparation). Table 
[2] gives the aperture magnitude measurements from the 
PAIRITEL observations. In this table and in all plots 
and modeling, we exclude exposures in which the CCD 
response was nonlinear, as well as -H-band observations 
during periods of poor transmission. 

We determined the position of the afterglow to be 
(a, 8) = (05 h 46 m 24.47 s ± 0.16", +64°12'38.60"± 0.18") 
(J2000) from a stacked J-band image covering the first 

30 min after the trigger. The quoted uncertainties are 
la, dominated by the astrometric mapping error from 
our stacked image to the catalog positions of 90 2MASS 
sources. 

2.5. NMSU 1 m Telescope 

Optical observations in the Johnson-Cousins UBVRI 
filters were obtained using the New Mexico State Univer- 
sity robotic 1 m telescope located at Apache Point Ob- 
servatory. Because the telescope happened to be pointed 
relatively near the burst location, the first observations 
were started only 47 s after the burst trigger, and only 

31 s after the alert. The telescope took five sequences 
of observations in the order /, R, V, B, U with exposure 
times of 10, 10, 20, 40, and 60 s, respectively, and then 
took another five sequences of observations with expo- 
sure times of 20, 20, 40, 80, and 120 s. With overhead, 
these ten sequences took about an hour to complete. The 
afterglow was detected in all of the observations, with 
random errors ranging from better than 0.01 mag in the 
first series to 0.1-0.2 mag in the last series. 

The afterglow brightness was measured using aperture 
photometry with an aperture of 3" radius; several refer- 
ence stars in the field were also measured. Calibration 
was achieved via observations of these reference stars, 

12 An animated image showing the variable trans- 
missio n and fading afterglow can be view ed at 
http://lyra.berkeley.edu/~jbloom/grb061126a-hband.mpg . 
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Fig. 1— False-color PAIRITEL finding chart (260" x 260") of the afterglow of GRB 061226 (left). The 2MASS catalog 
stars used for the photometric calibration are denoted with boxes at the catalog positions and labeled with J-band (2MASS) 
magnitudes. At right, IR images show the fading of the afterglow from 58 s to 1 d after the GRB; the images are progressively 
deeper at later epochs. 



TABLE 4 

KAIT and Lick 1 m Observations of GRB 061126 



TABLE 5 
UVOT Observations of GRB 061126 



Ai a 


Filter 


Exp. Time 


Magnitude 13 


Flux b 


00 




w 






366.0 


clear 


20.0 


16.29 ±0.12 


946.2 ± 99.0 


409.5 


V 


45.0 


17.25 ±0.12 


464.5 ± 48.6 


465.5 


1 


45.0 


16.24 ±0.11 


776.2 ± 74.8 


3336.0 


R 


300.0 


18.86 ± 0.05 


88.7 ±4.0 



Note. — i?-band points are from the Lick 1 m 
(Nickel) telescope; all others are from KAIT. Only the 
first exposure in each filter is given; additional data are 
available in the online supplement, or by request. 
a Exposure mid-time, measured from the Swift trigger 
(UT 08:47:56. 4). b Observed value; not corrected for 
Galactic extinction. 

along with UBVRI standard stars, on 22 December 2006. 
A standard photometric solution was derived for this 
night, yielding calibrated magnitudes for the reference 
stars. On the calibration night, the standard-star solu- 
tion yielded rms deviations of about 0.03 mag in each 
bandpass, so the calibration zeropoints are accurate only 
to this level. The final photometry is given in Table El 

2.6. KAIT and the Lick Nickel Telescope 

The 0.76 m robotic Katzman Autom atic Imaging Tele- 
scope ( KAIT; iFilippenko et aTll2001h and its GRB alert 
system (|Li et al.l l2003bl ) responded to the GCN notice 
with the position of GRB 061126 at t = 16 s after the 
trigger, and attempted to execute a pre-arranged obser- 
vation sequence. Unfortunately, the weather conditions 
were poor, so KAIT did not acquire a useful image un- 
til t = 305 s. Also, the telescope did not have a suc- 
cessful focusing procedure before the GRB observations 
(again due to bad weather) , so the images were not fully 
in focus. Nevertheless, a sequence of V, I, and unfil- 
tered observations was made, and the GRB afterglow 
was detected in most images. A successful focusing pro- 
cedure was executed during the middle of the GRB ob- 
servations, and KAIT followed the GRB until humidity 
forced the system to shut down at t ~ 1.8 hr. The images 



At a 


Filter 


Exp. Time 


Magnitude 13 


Flux b 


00 




w 


(MJy) 


2152.0 


U 


967.5 


18.64 ±0.11 


67.1 ±6.6 


2114.0 


B 


809.5 


19.05 ±0.12 


99.2 ±10.1 


2178.0 


V 


809.5 


18.43 ±0.14 


156.5 ± 19.1 


8862.0 


V 


902.0 


20.51 ±0.32 


23.0 ±5.8 


7700.3 


u 


196.6 


19.75 ±0.19 


24.2 ±3.8 


15119.5 


u 


295.1 


20.33 ± 0.38 


14.2 ±4.2 


2807.3 


B 


196.6 


19.52 ±0.12 


64.4 ±7.0 


7905.3 


B 


196.6 


20.61 ±0.26 


23.5 ±5.0 


3421.3 


V 


196.6 


19.07 ±0.23 


87.2 ±16.5 



Note. — UVOT data points were not used in our 
light-curve models. U, B, and V filter measurements 
are our own re-reductions. UV and white filter pho- 
tometr y are from GCN Report 16.2 iSbarufatt i et al.l 
2006b); they are not repeated here, but are also avail- 
able in our online supplement. 

a Exposure mid-time, measured from the Swift trig- 
ger (UT 08:47:56.4) . b Observed value; not corrected for 
Galactic extinction. 



were automatically processed with the proper dark cur- 
rent, bias, and flat fields before measuring photometry 
and calibrating relative to ten SPSS stars (|Coo]| [20061 
lAdelman-McCarthv et~aT1l2006l ). The reference magni- 
tudes of these stars were converted to VRI using the 
transformation equations of Lupton (2005) 13 . Unfiltered 
observations were calibrated using the i?-band magni- 
tudes. 

A sequence of 300 s i?-band images was also manually 
obtained at the Lick Observatory Nickel 1 m telescope 
from t w 1.0 hr to 1.7 hr. The observations were again 
terminated prematurely due to the weather conditions. 
The images were manually reduced with the proper cal- 
ibration files (bias, dark current, and flat-field images), 
and calibrated relative to our SDSS reference stars. The 
photometry is given in Tabled 

13 http: / /www.sdss.org/dr4/algorithms/sdssUBVRITransform.html#:Luptonl 
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2.7. Swift UVOT 

Swift began follow-up observations of GRB 061126 at 
1605 s after the burst trigger, its slew having been de- 
layed as a result of the Earth-limb constraint. Despite 
this time delay the afterglow was still detected, albeit 
marginally, in all of the UVOT filters except for UVW2. 

We acquired the UVOT imaging data from the NASA 
archive 14 . Unfortunately the afterglow had already be- 
come quite faint by the beginning of the observations, 
and is not well detected except in stacked exposures. 
To calculate the most accurate photometry possible, we 
therefore bin observations obtained at similar epochs, 
and perform aperture photometry using the optimal 
apertu re size as given in the prescription of iLi et all 
(|2006h for the V, B, and U data. For the UV filters, we 
use the photometry reported by the U VOT team in the 
most r ecent GCN Report for this burst (fSbarufatt i et al.l 
2006b]). Our photometry is given in Table [5] 

Despite our small-aperture analysis, the UBV photom- 
etry has large uncertainties, and some points are not for- 
mally consistent with effectively simultaneous ground- 
based observations in the same filters. This is not neces- 
sarily surprising; due to the extremely faint nature of the 
afterglow by the time that Swift completed its slew, the 
actual uncertainty on these measurements may be signifi- 
cantly larger than the nominal photometric error. We do 
not use these points in our modeling, but we do include 
them in our plots. The UVW1 and UVM2 points have 
even larger uncertainties (0.4-0.8 mag), and no ground- 
based calibration is available; moreover, Galactic extinc- 
tion (which is significant in this direction: Ajjvm2 ~ 1.8 
mag) is increasingly uncertain toward these wavelengths. 
Therefore, we exclude these points from the formal fits as 
well, and restrict our modeling to the much more precise 
ground-based photometry. However, we do include the 
UV points in our SED plots for comparison. 

2.8. GCN Circulars 

For comparative purposes, our plots also include points 
from the GCN Circulars 15 , but we do not actually use 
these points in our fitting. Most early-ti me data were 
calib rated against the USNO Bl.O survey (jMonet et al.l 
|2003[ ) , which does not contain very accurate photometry. 
Some data points were also calibrated against a prelimi- 
nary release of an SDSS pre-burst observation of this field 
(ICool et al.ll2006T) which was later found to be incorrectly 
calibrated (indeed, our own use of these observations for 
our preliminary calibrations exposed the problem and 
motivated the re-releas e of the SD SS calibration used for 
the KAIT reductions, ICooll 120061 ) The photometry is 
summarized in Table [5] 

2.9. Keck Host Imaging and Spectroscopy 

The galaxy hosting GRB 061126 was observed on 18 
January 2006 wit h the Low Res olution Imaging Spec- 
trometer (LRIS; lOke et al.l I1995Q on the Keck I 10 m 
telescope. A total of 3 x 300 s of imaging was obtained 
simultaneously in R c and g' filters, followed by two long- 
slit spectroscopic exposures of 20 min each over the host 
location identified in our images (Figure [2|). 



Imaging data were reduced using standard techniques. 
In the final, reduced images we identify an extended 
source consistent with the GRB afterglow position, ap- 
parent in both our g' and R imaging. This source, which 
we identify as the GRB host, has significant substruc- 
ture, with two clearly distinct elongated knots, each of 
which may have further substructure. The afterglow is 
offset by < 1" from the brightest knot. The magni- 
tude of the whole complex, using a large (2.55" radius) 
aperture fixed at the position (a, 5) = 05' l 46 m 24.53 s , 
+64°12'39.31") and performing photometric calibration 
using the subset of our SDSS stars that are unsatu- 
rated in the Keck images, is g = 25.13 ± 0.16 and 
R = 24.10 ± 0.11 mag. After correcting for Galactic 
extinction, the g — R color is 0.88 ± 0.19 mag. 

Using the long slit of width 1.0", the 600 line mm -1 
grating blazed at 7500 A, and the 300 line mm -1 grism 
blazed at 5000 A, we obtained spectra covering the wave- 
length range 6680-9266 A and 3200-7649 A with the red 
and blue cameras, respectively. The spectroscopic data 
were processed with an IDL package 16 customized for 
LRIS long-slit reductions developed by J. Hennawi and 
J. X. Prochaska. The observed PA of -69.35° is signifi- 
cantly different from the parallactic angle, but as we are 
only seeki ng a line identifica tion the effects of differential 
slit losses (|Filippenko 1119821 ) are not very significant. 

The two-dimensional reduced spectra show a faint blue 
continuum and a sole emission feature at a vaccuum 
wavelength of ~ 8050 A (Figure [3]). From the one- 
dimensional spectrum it is clear that the emission feature 
is slightly resolved; it has a full width at half-maximum 
intensity (FWHM) of roughly 10 A, while the spectral 
resolution in that wavelength range, as determined by 
the FWHM of arc-line profiles, is 5.8 A. As we detect 
only one emission feature, it is impossible to definitively 
report the redshift of the GRB host. However, the width 
of the line and the presence of obvious continuum blue- 
ward of this feature suggest that it likely corresponds 
to the [O II] AA 3727.11, 3729.86 A doublet. Using our 
best-fit observed line centroid of 8049.0 ± 2.0 A and as- 
suming an intrinsic doublet centroid of 3728.5 ± 0.5 A, 
this would imply that the GRB host lies at a redshift of 
z = 1.1588±0.0006. The width of the emission feature is 
consistent with an [O II] doublet at this redshift (3.19 A 
separation at z = 1.159) when convolved with the instru- 
mental resolution. Furthermore, if the emission is from 
[O II] we would not expect to see any other spectral lines; 
the most common, redder lines would fall redward of our 
red spectrum, and the only strong emission line blueward 
of [O II] is Lya which lies below our spectral coverage. 

Adopting the 8050 A emission feature as the [O II] 
doublet, we can estimate the star-formation rate (SFR) 
of the host galaxy based on the line flux. We measure the 
total flux from the extracted, one-dimensional spectrum 
in the region 8042.0-8056.0 A. We measure the contin- 
uum by computing the median flux in two regions free 
of sky lines (regions "C" in Figure [3]). The total flux 
in these lines is (1.6 ± 0.2) x 10~ 17 erg s^ 1 cm -2 . Us- 
ing the rel ation betw e en [0 II] luminosity and SFR de- 
scribed bv lKennicuttl (|1998fl . we find that the GRB host 
galaxy is undergoing star formation at a rate of 1.6 ± 0.2 



14 http:/ /heasarc. gsfc.nasa.gov/docs/swift/archive/ 

15 http://gcn.gsfc.nasa.gov/gcn3_archive.html . 



16 Sp ecifically, the Long slit codes now bundled within XIDL, 
http://www.ucolick.org/~xavier/IDL/indcx.htnd] . 
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Fig. 2. — Keck I (LRIS) late-time finding chart. The image is 25" on a side and represents the stacked g' and R bands. 
The placement of the l"-wide slit is shown in green. The white circle is the 2a position of the afterglow as measured from 
PAIRITEL images. The probable host galaxy is visible at this location. The object is blue and extended, appearing to have a 
complicated morphology. Further down the slit 21.1" to the southeast is a relatively bright galaxy with z = 0.6225 ± 0.0004, 
based on emission from [O III] AA4960.2, 5008.2, H/3, and [O II] AA3727.11, 3729.86. (This second galaxy is not associated with 
the GRB.) 



M© yr -1 (assuming H — 71 km s _1 Mpc -1 , f^M = 0.3, 
= 0.7, as is done throughout this paper). The intrin- 
sic uncertainty in the calibration converting JO II] line 
lumin osity into SFR is approximately 30% (jKennicuttl 
1998); this uncertainty is not included in our quoted er- 
ror. Though our measured SFR must be considered a 
lower limit because we have not accounted for dust ex- 
tinction, it is interesting to note that the rate of star for- 
mation in the host galaxy of GRB 061126 is comparable 
to the SFRs found in other GRB hosts, measured using 
dust-corrected UV fluxe s as well as [O II] luminosities 
(|Christensen et al.ll2004) . 

3. DATA ANALYSIS AND MODELING 
3.1. Prompt Emission 

The prompt emission from GRB 061126 displays a 
complex, multi-peaked profile dominated by two large 
pulses (Figure [4}. After a gradually increasing compo- 
nent starting at 7 s before the trigger and a small pre- 
cursor spike at t ~ 1 s after the trigger, the first and 
largest pulse ("Pulse A") begins at t w 3 s, and fades 
to about twice the background level by t « 15 s. A 
second pulse ("Pulse B") begins at i « 19 s, lasting 
until i « 25 s. There is short-timescale microstruc- 
ture in both of these p ulses. The full burst has a Tgo 
(|Kouveliotou et al.|[l993j) of 26 8 ± 8 s in the full Swift 
15-350 keV band (|Butler et al. l l2007f) . 

Using the combined Swift BAT and R HESSI spectrum, 
we fit a Band model (|Band et al.lll993T ) over several dif- 
ferent time regions: The entire burst, the first pulse, and 
the second pulse. Results are summarized in Table[7J For 
the full burst, we measure an average peak energy £^ j0 f)s 



of 620 keV, though there probably is hard-to-soft spec- 
tral evolution, as evidenced by the significantly different 
values of E Pj0 b s during the two pulses. The total fluence 
over the full spectral range is (3.0±0.4) x 10 -5 erg cm -2 , 
which at the putative host redshift of 1.16 corresponds to 
an isotropic release of energy of Ei SO = (1.06±0.14) x 10 53 
erg over a 1-1000 keV host-frame bandpass, assuming 
our standard cosmology. 

Given this val ue of Ei SO , the Amati relation 
(jAmati et al.ll2002T ) predicts E Pi0bs = 130 x 10 0±0 - 3 keV, 
which is 2cr from the measured value. Given Ej K „ and 
Ep^obs, the Ghirlanda relation (Ghirla nda et al.l 1 20041 ) 
predicts a jet break at t = 10 141±0 - 27 d (13-49 d). The 
beaming-corrected energy release is 10 51 - 6±0 - 2 erg, for a 
jet opening angle of K) 1 - 2 ^ 010 degrees (13-20°). 

RAPTOR detected GRB 061126 contemporaneously 
with the BAT emission. Fortuitously, the first unfiltered 
exposure (which took place 20.87-25.87 s after the BAT 
trigger) matches quite well the second pulse of the GRB 
(peak time of 22.5 s and a Tqq of about 5 s) — see the 
first optical point of Figure [4] 

Comparing this first, contemporaneous data point to 
later exposures in which the prompt emission has faded, 
the early-time RAPTOR data are seen to fade as a sim- 
ple power law [F v oc i -Q ), with the decay index a « 1.5 
(using the GRB trigger time as to). There is no evidence 
for an additional prompt flux component based on an ex- 
trapolation backward from later measurements. Consis- 
tent w ith this observati on, if we extrapolate the best-fit 
Band (|Band et al.| [l993) model of the second GRB pulse 
into the optical (Figure[5]), the predicted optical flux is 
only 250 //Jy, significantly less than what is observed. 
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Fig. 3. — Spectrum of the host galaxy of GRB 061126 in the vicinity of the sole emission line, interpreted as the [O II] 
doublet. The data were taken with LRIS on Keck I using the 1" slit and the 600/7500 grism. Top: Background subtracted, 
two-dimensional spectrogram showing the weighted mean of two 1200 s exposures. Pixels contaminated by cosmic rays in one 
exposure are excluded from the mean. Bottom: One-dimensional, coadded spectrum. The profile fit from a bright, nearby source 
was used to extract a spectrum at the known location of the GRB host galaxy. The dashed line represents the la uncertainty 
at each pixel. The width of the emission feature is comparable to the spacing of the [O II] doublet at a redshift of 1.16. The 
doublet spacing is approximately equal to the instrument's resolution element at 8050 A, so we expect the doublet to be barely 
resolved. We measure the galaxy's continuum by computing the median flux in two regions free of night-sky lines (regions "C"). 
Using this continuum value, and measuring the signal between 8042.0 and 8056.0 A, we find that the total flux in this emission 
feature is (1.6 ± 0.2) x 10~ 17 erg s" 1 cm" 2 . 
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Fig. 4. — BAT 15-350 keV light curve of GRB 061126, showing the profile of the main burst (dominated by two primary pulses, 
labeled "A" and "B") as well as the rapid response from the robotic telescopes RAPTOR, the NMSU 1.0 m, and PAIRITEL. 
At this early time the light curves are well fit by a simple power law {F v oc t~ 1,5 ) with to simply set to the trigger time. There 
is no evidence for a rising component or any correlation of the optical emission with the prompt emission. 




Fig. 5. — High-energy SED during the second main pulse of the prompt emis sion ("Pulse B"), with the contemporaneous 
RAPTOR measurement superimposed. An extrapolation of our best-fit Band ( |Band et al.l 119931 ) model (Table [7]) into the 
optical is shown to underpredict the optical flux by several orders of magnitude, suggesting that even at this very early time 
the optical afterglow is already present at 20 s post-trigger, dominating the early-time flux at long wavelengths. A temporal 
extrapolation of the BAT light curve to 60 s similarly underpredicts our multi-color data at that time. 
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Finally, if we extrapolate the Band model of the prompt 
emission in time to 60 s (assuming continued hard-to-soft 
spectral evolution and fading), it both falls far short of 
our multi-color data at that time and also has a different 
spectral slope (open circles and dashed line in Figure E|) . 
For this burst, there is no clearly observable association 
between the prompt emission and the long-wavelength 
afterglow, even as early as 20 s after the burst. This 
is not inconsistent with earlier reports of a link (e.g . , 
iBlake et alJl200H iVestrand et aT1l200& lYost et al.l[200% : 
more likely, as our SED shows, for this burst the prompt 
component is simply dominated by an extremely bright 
early-time afterglow. 

3.2. Optical Light Curve 

While the very early optical light curve appears to fol- 
low a power-law behavior, the light curve enters a more 
complex phase within a few minutes. A brief visual in- 
spection of the overall afterglow light curve (Figures [6] 
and [7]) shows several interesting features of this burst. 
The early decay that began in our earliest data continues 
with the same trend for several decades in time, fading 
roughly as a power law with an average decay index a 
of about 1.5. A simple power law is clearly a poor fit, 
however, due to the presence of of a prominent "bump" 
feature at around 120 s. 

Rapid early decay gives way to a shallower-decaying 
(a < 1) component starting at about 10 3 s that dom- 
inates for the rest of our observations. Reports from 
the GCNs indicate an apparent stee pening at later 
times, though with d i fferent auth ors (|Rol et ail 120061 : 
iPozanenko et all 120061 : iMisral l2006f) disputing the value 
of the late-time decay index. 

In our modeling, we construct the light curve as a sum 
of several broken power-law components, each of gen- 
erally different color as well as different decay indices. 
Details of the model follow. 

The broken power -law model is that of 
iBeuermann et all (|1999f ). Mathematically, the fit- 
ted curve (as a function of the filter, represented by a 
central frequency v) has the form 

The flux offset F V j generally depends on the filter (rep- 
resented by the subscript v) as well as on the component 
j. The temporal decay index a depends on the compo- 
nent and on whether it is before the break (subscript b) 
or after (subscript a), but does not depend on the fre- 
quency. 

In our modeling, we use three total summed compo- 
nents, though two of them are tightly linked: a "fast" 
component, with no break (a^ = oiF,b), a "bump" com- 
ponent, fixed to have the same color as the fast compo- 
nent {F v> b oc F Vi f), and a "slow" component, also with 
no break (as. a = czs,b)- 

The constraints on the slow component depend on the 
specific model tested. In all cases, the decay index (ccg) 
is free to vary. We do not directly constrain the spectral 
index (defined here using the convention F v oc v~fj) or 
require a power-law SED. However, in some models we 



do constrain the change in the spectral properties. In our 
different models, these color constraints are as follows: 

• Unconstrained, allowing the linear flux factor F„ t s 
to take arbitrary values for each filter. 

• Constrained to allow no spectral evolution: F Uj s = 
CFv^, where C is a fitted constant. 

• Constrained in such a way as to permit the spec- 
tral index of the slow component to differ from 
that of the fast component by A/3p-s by fixing 

F v ,s = CF v ,F(y/^R-band)~ APF - s ■ Note that this 
difference in the spectral index can be fit directly 
regardless of Galactic extinction (which is signifi- 
cant) or host-galaxy extinction as long as the ex- 
tinction can be assumed to be a constant in time. 

We also attempt two other, somewhat more specialized 
models: 

• A fourth, physically motivated model in which the 
slow component is treated as an adiabatic forward 
shock whose flux peaks at some time during our 
multi-color observations. In this case the slow com- 
ponent is broken, where the rising and decay in- 
dices (as,b and ocs,a-> respectively) and the change 
in spectral index across the break (A/3g) ar e fixed 
using the constraints from the X-ray spectrum dis- 
cussed in i j4.ll 

• Solely for the purposes of estimating the "average" 
early-time optical decay index, we also fit a model 
with no bump component at all. (This is a very 
poor fit, and because most of the IR observations 
were performed during the bump any photometric 
SED generated would be highly biased. Therefore 
we do not use this model for any other purpose.) 

Our data only trace the light curve with no gaps in 
coverage until ~ 7000 s, and we have no color informa- 
tion past 4000 seconds aside from a marginal J-band 
detection and H and K s limits from the second night. 
Later-time measurements are present in the GCN Circu- 
lars, but different authors report different and somewhat 
contradictory behavior, suggesting either a problem with 
some of the public data or complex behavior of the light 
curve at late times. Consequently, we do not include 
any optical points after 10 4 s in our fits, anticipating 
that the late-time optical evolution will be discussed in 
greater detail in upcoming work by Mundell and collab- 
orators. Instead, we focus on the properties of the early 
and intermediate afterglow. 

After fitting, we correct for Galactic extinction of 
Eb-v = 0-182 mag (using the NED extinction calcu- 
lator 17 , ISchlegel et al. 1998) and fit a simple power law 
to estimate the observed spectral index for different com- 
ponents. (This neglects the possibility of host extinction, 
but as we will show in H4.41 if host extinction is present 
it does not cause significant deviation from a power law.) 

17 The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space 
Administration. 
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Fig. 6. — _R-band and X-ray light curves of the afterglow of GRB 061126 showing the basic features of the early-to-late optical 
and X-ray afterglow light curves. Unfiltered data are also included, offset by 0.13 mag to match the _R-band calibration. Optical 
data are shown as filled (used in modeling) or unfilled (not used) symbols; X-ray data are shown as error bars with no central 
symbol. There is a rapid decay with a bump at early times, transitioning to a significantly slower decay that probably breaks 
at late times. Due to a delayed slew the X-ray afterglow was not observed until 1600 s, but from then until it faded below the 
detection threshold at ~ 9 d it decays as an unbroken power law. The 1 keV normalized flux is plotted. 




Fig. 7. — Broadband light curves of our early-time, multi-color photometry of the afterglow. The light curves are fitted with 
a three-component broken power-law model, assuming that the third component (with the slowest decay) has a spectral index 
that differs from that of the early component by an amount A/3. 
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Fig. 9. — Evidence for color change across the transition from the rapidly decaying component to the slowly decaying component. 
We fit a power law to simultaneous J/K B PAIR1TEL exposures and optical exposures (in U, B, V, R, or /; color-coded 
appropriately) from the NMSU 1.0 m telescope. The solid curves are not direct fits to these data, but represent the spectral 
index that would be observed at each time if fit to K s — U photometry based on two of our models. The solid line is for a 
model where the late-time component of the afterglow is modeled as a simple power law; the dashed line represents a model of 
a forward shock undergoing a minimum-energy break at approximately the transition time. 



Despite the complexity of the models, no fit is observed 
to give a value of x 2 P er degree of freedom (dof) w 1. 
The most successful model, using unbroken power laws 
and arbitrary filter-dependent color change, gives x 2 /dof 
w 6. This is not surprising; modulations in afterglow flux 
have been observe d in many previous cases (for example, 
iLipkin et al.l 120041 ) . Accordingly, we add in quadrature 
an extra uncertainty of 0.08 mag to all the photometric 
measurements. With this adjustment, our different mod- 
els produce fit results as given in Table [Sj There is strong 
evidence for color change — a fit allowing no change in 
the relative fluxes of different components has a value 
of x 2 /dof = 261.9/193, while an equivalent fit allowing 
the early fast-fading and later slow-fading components to 
have different spectral indexes (Figured]) gives x 2 /dof = 
231.3/192. (The /-test significance of this improvement 
is 10~ 6 .) Under this model, (3 shifts toward the blue 
across the transition by Af3p-s — —0.383 ± 0.075. 

3.3. Optical Color Evolution 

The above results indicate color evolution across the 
steep-to-shallow transition. To confirm this behavior and 
ensure this is not an artifact of the modeling, we con- 
struct truly contemporaneous SEDs using simultaneous 
PAIRITEL and NMSU 1.0 m observations by mosaicing 
only those individual PAIRITEL exposures taken during 
the NMSU exposure time range. In cases where the ex- 
posure time is much less than the time since the burst, 
we increase the effective PAIRITEL exposure time by 
also including exposures shortly before the beginning or 
after the end of the optical exposures (otherwise the S /N 
for the late-time PAIRITEL measurements is quite low) . 



After correcting for Galactic extinction, we fit the three 
J/ .ftf s /optical data points with a basic power law, and de- 
termine the best-fit value and uncertainty for the spectral 
index /3. (We omit _ff-band measurements from this fit, 
since the variable transmission introduces a significantly 
larger scatter in that band compared to J and K s .) The 
results are plotted in Figure [9j 

The early-time colors from the first NMSU filter cycle 
are consistent with a spectral index of (3 — 1.2 ±0.1, but 
starting at around 500 s the colors shift notably blueward 
and at later times the index is typically (3 — 0.95 ± 0.10. 
This provides model-independent support of our fit con- 
clusion that the afterglow has undergone a color change. 

3.4. X-ray Light Curve and Spectrum 

The XRT light curve (Figure [6]) fades as a purely un- 
broken power law with a decay index ax = 1-31 ± 0.01 
over the entire span of the Swift observations, from ~2 ks 
out to nearly 10 d. There is some suggestion that the de- 
cay rate is slightly faster during the first orbit, but the 
unbroken fit is still good (x 2 /dof = 668.2/550) and joins 
with an unbroken extrapolation of the fading BAT light 
curve. 

After removing the effect of neutral hydrogen 18 absorp- 
tion [both Galactic, for which we estimate A^h, Galactic = 
0.103 x 10 22 cm" 2 from lDickev fc Lockmanl (|1990f l. and 
at the host redshift, for which we calculate a best-fit value 

18 More accurately, the X-ray absorption is dominated by light 
metals, rather than hydrogen. However, we will use the standard 
terminology and refer to the absorption column in terms of the 
hydrogen density. 
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of A?H,ho8t = (1-1 ±0.3) x 10 22 cm -2 ], the X-ray spectrum 
is a good fit (x 2 /dof = 213.9/238) to a simple power law, 
with a spectral index (3x = 1-00 ± 0.07. There is no 
evidence for spectral evolution durin g the observations: 
using the X-ray hardness ratio (see iButler fc Kocevs"kH 
l2007f ): we constrain the change in the X-ray spectral 
slope to be Aflx < 0.4 (90% confidence). Our analy- 
sis is consistent with that from the most r ecent GCN 
report for this event (|Sbarufatti et al.ll2006bl ). 19 

4. DISCUSSION 
4.1. Constraints on Synchrotron Model Parameters 

Due to the complex behavior of the optical light curve 
and the possibility of host extinction, it is difficult to ap- 
ply any constraints on the physical parameters of this 
burst (p, cb, £_b, etc.) from the optical observations 
alone. However, in principle the degeneracy can be bro- 
ken using the XRT data. 

The fits to our data unambiguously show that the op- 
tical flux and X-ray flux are decaying at different rates 
during the period when we have contemporaneous data 
in both bands (1800-6000 s): a opt = 0.75 ± 0.05 and 
ax = 1.31 ± 0.01 . Assu ming a normal synchrotron spec- 
trum (Sari et al. 1998), this allows only two possibili- 
ties for the ordering of the spectral breaks: v c < v opt < 
v m < fx (fast cooling), or v m < v opt < v m < v x (slow 
cooling). In either case, the X-ray spectral index is pre- 
dicted to be (3 — p/2. The observed spectral index is 
= 1.00 ± 0.07, so we can confidently calculate an elec- 
tron index of p — 2.00 ± 0.14, the minimum value for a 
distribution that is unbroken at high energies. 

However, this conclusion is problematic in light of the 
observed X-ray light curve. Independent of environment, 
the decay index of emission from an expanding adiabatic 
forward shock in the regime v c , v m < v should obey a = 
3/3/2 — 1/2, which for our spectral data would imply a = 
1.00 ± 0.10, disagreeing by about 3tr from the measured 
value. 

The discrepancy is significant and indicates that at 
least one standard assumption does not hold. We con- 
sider several possibilities: 

1. Our conclusion of u c ,v m < vx is incorrect. But, 
this is unlikely: no other ordering can produce op- 
tical and X-ray light curves that decay at different 
rates. Furthermore, the X-ray light curve remains 
unbroken out to very late times ( > 10 6 s). Such a 
late cooling break in the X-ray band would imply 
an extremely low density (-C 1 cm -3 ), inconsistent 
with an intragalactic environment and the high Nh 
column density observed in the X-ray data. (A col- 
umn depth of 1 kpc predicts n > 4 cm -3 if the gas 
is uniformly distributed throughout the galaxy.) 

2. The afterglow evolution is dominated by radiative 
losses and has not yet transitioned to an adiabatic 
phase, even out to ~ 9 d. This case would pre- 
dict that a = 12/3/7 - 2/7 = 1.42 ± 0.12, within 
la of the observed value. However, this requires 

19 Note, however, that the GCN report quotes the photon spec- 
tral index, instead of the flux spectral index, and that since they 
fit neutral hydrogen absorption only at z = instead of the host 
redshift, their Nh excess is significantly less than our TVjj host- 



that the synchrotron spectrum still be fast-cooling, 
which demands that the optical band be above the 
synchrotron critical frequency [y c < v opt < u m ) to 
produce a fading optical light curve in a constant- 
density medium. The specific, predicted value of 
a = 4/7 = 0.57 is not statistically consistent with 
our observations of the optical decay index in this 
region, but it is possible that the measurement 
may be confused somewhat by the presence of soft 
breaks or other systematic effects. 

3. The X-ray light curve is affected by synchrotron 
self-Compton (SSC) losses, which would steepen 
the decay slightly relative to that pre dicted by a no- 
Compton model (|Sari fc Esinll200lb . The amount 
of predicted steepening is, unfortunately, relatively 
low: Aa « 0.5(p - 2)/(4 - p), which for the X- 
ray inferred value of p = 2/3 = 2.0 ± 0.14 gives 
a constraint of the SSC-corrected decay index of 
1.0 ± 0.14, still off by 2a from the observed value 
of 1.31 ±0.01. 

4. The X-ray emission is not due to a forward shock, 
or not due to synchrotron radiation. 

Some combination of factors may also be at work. In 
any event, the specific resolution to this problem is not 
critical at this stage, and we will proceed assuming p = 2, 
which is independent of any assumptions about the burst 
environment and requires only that the source of the X- 
ray light be synchrotron radiation, where the X-ray band 
is above the cooling and minimum-energy breaks. 

In addition to identifying the spectral regimes, we also 
seek to identify the burst environment (constant-density, 
wind-stratified, or a spreading jet). A spreading jet has 
a rapidly decaying light curve and can be immediately 
ruled out. A wind can also be ruled out, since it would 
predict that the optical flux should decay faster than 
the X-ray flux, yet we observe the opposite. However, 
a constant-density (interstellar medium) environment is 
consistent with all of the qualitative features of the light 
curve. 

4.2. The Very Early Afterglow Decay — A Reverse 
Shock? 

For any assumption about the X-rays, the behavior 
of the fast-decaying optical component cannot be repro- 
duced as originating from a forward shock. The early op- 
tical decay index is Qf ss 1.5 if we fit a single component 
(no bump), though in our fits with a bump component 
the decay rate is steeper. (However, if we interpret the 
bump as a density variation or other modulation of a sin- 
gle power law, rather than a truly separate component, 
the value of cxf s=s 1.5 is most appropriate. For our dis- 
cussion of the interpretation, we use the constraint a > 
1.5, which is certainly true independent of the number of 
components used in our model.) 

This decay is steeper than the ax = 1.31 decay slope 
observed in the X-rays at late times. In no regime of the 
forward-shock interstellar-medium model can the longer- 
wavelength emissio n decay faster th an the X-rays at the 
same or later times (jSari et al.lll998[ ). 20 Furthermore, the 

20 It is possible in a wind-stratified medium, but even then, the 
observed value of cip is too steep for p = 2. 
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subsequent passage of v m or v c through the optical bands 
should produce breaks in the optical light curves which 
can only accelerate the temporal decay, not decelerate it. 

The lack of flux excess in the first RAPTOR mea- 
surement, and the extrapolation of the Band model 
into the visible frequency range, argue against associa- 
tion with this emission to a prompt optical component 
from the internal shocks that (presumably) produced the 
high-energy emission. Therefore, the most likely candi- 
date to be responsible for this early optical emission is 
an external reverse shock which is produced when the 
forward shock begins to interact with the surrounding 
medium. Such an explanation has been previously used 
by several authors to explain observed peculiarities in 
the ea rly-time light curves of a number of GRB after- 
glows dSari fc Pirarifl999 t lKobavashi h Zhand l2003l; IWeil 
12001 IShao fc paill2005f) . 

iKobavashil (|2000f ) derived the temporal behavior of 
reverse-shock light curves for two different regimes: the 
thick shell (corresponding to cases in which the reverse 
shock is able to become relativistic as it crosses back 
through the ejecta) and the thin shell (for cases where 
the shock always remains nonrelativistic). The thick shell 
case is associated with relatively long bursts, and in par- 
ticular predicts that the peak of the reverse-shock emis- 
sion will happen at about the same time as the burst it- 
self, whereas the thin-shell case predicts an optical peak 
after the GRB. As the RAPTOR data indicate an op- 
tical light curve that is already rapidly declining during 
the prompt emission, we apply the thick-shell case here. 
(However, the thin-shell case gives identical numerical 
predictions for the region of interest for p = 2.) 

Before the shock crossing time T, the reverse-shock 
light curve is expected to rise as t 1 / 2 , which as already 
discussed is not observed. The light curve peaks at T 
and then decays: for frequencies v m < v < v c the flux 
decays as £-( 73 p+ 21 )/ 96 ; while for v < v m it decays as 
£ 17 / 36 . (For v c < v the light curve shuts off and rapidly 
dives to zero, and for fast cooling the light curve either 
declines asa= 17/36 or is shut off after T .) 

Observationally, the early-time average decay index is 
otF,av ~ 1.52 ±0.02 if we fit only a single component with 
no bump, though because of the bump feature the exact 
decay rate clearly varies; given this complexity, a rigorous 
statistical comparison with the simplistic predictions is 
not possible. However, we can unambiguously rule out 
all cases except for the slow-cooling case in which v m < 
v < Vc, where, if we use the value of p from the late-time 
X-ray spectrum, we predict a — 1.74 ± 0.11, which is 
generally consistent with our observations. 

Our early-time multi-color photometry provides us 
with a rare opportunity to test the spectral predictions 
of reverse-shock models, in addition to the time-domain 
predictions. The synchrotron spectrum of a reverse shock 
is predicted to be identical in its general form to that 
from a forward shock, with the exception that after the 
reverse shock crosses the shell (at time T), there is no 
emission from frequencies above some cutoff frequency 
is cu t- The spectral power- law indices are the same (for 
the same value of p) , although the break frequencies v m 
and v c are generally not. 

For the reverse shock, the only regime consistent with 
the data is v m < v < v c . The prediction for the spectral 



index in this region (from|Sari_e£al. 1998) is v Kl a . 
Again assuming p = 2.00 ±0.14 from the late-time X-ray 
observations, we predict that the spectral index should 
be p = 0.50 ± 0.07. This is significantly less than the 
observed value (after correcting for Galactic extinction) 
of Pf w 1.25. This very large discrepancy may be due 
to host-frame extinction. However, as we will discuss 
in £|4.4[ the lack of observable curvature in the SEDs of 
either the fast-decaying or the slow-decaying components 
impose strong constraints on the amount and nature of 
any extinction. 

4.3. Transition to Forward Shock 

The fast-decaying early-time optical component is sub- 
sumed by the shallower late-time component at around 
10 3 s. The natural inclination is to assume that the 
reverse-shock emission has been replaced by that from 
a forward shock. 

As already discussed, the different decay rates of the X- 
ray and optical light curves rule out the possibility that 
they are in the same synchrotron regime, so we expect 
the optical band to be below a spectral break. For adia- 
batic evolution, if the optical band is below the cooling 
break, we expect (assuming p — 2.0) an optical decay 
index of a op t = 0.75. This is consistent with the data, 
both in the empirical model where we assume a power 
law that fades throughout the observations (for which 
we calculate as = 0.75 ± 0.06) and the more realistic fit 
in which the slow component rises, experiences a (chro- 
matic) break, and then fades (x 2 /dof = 228.0/197 for 
this model, slightly better than the best model assuming 
an unbroken slow decay.) For radiative evolution and the 
case that the optical band is below the minimum-energy 
break (but above the cooling break and the critical fre- 
quency), a op t = 4/7, which is only consistent within 
2.5cr. 

If the slower decay component is due to an adiabatic 
forward shock and the optical band is in the range v m < 
v D pt < v c then the predicted intrinsic optical spectrum 
is F v oc ^(p -1 )/ 2 , which for p — 2.0 ± 0.14 suggests 
= 0.5 ± 0.07. If the forward shock is radiative, the 
prediction is the same but with no statistical uncertainty 
{(3 = 1/2). The observed spectral index of the possible 
forward-shock component is somewhat redder than this 
value, Ps = 0.92 ± 0.08 — or, using our physical model, 
the post-peak (3$ = 1-11 ± 0.09. Again, one may appeal 
to host-galaxy extinction to make up the difference. 

We will examine this possibility in greater detail in 
the next section. For now, however, we can sidestep the 
host extinction question completely if we look only at 
the change in spectral index between the early (reverse 
shock) and late (forward shock) components, as long as 
we can assume that the extinction is constant with time. 

As discussed earlier, the only reverse-shock model con- 
sistent with the observed early-time light curve requires 
that either v m < v opt < v c (adiabatic) or v c < is opt < v m 
(radiative), both of which predict the same value for the 
spectral index for p w 2. We reach the same conclusion 
with the later-time data, meaning that the reverse and 
forward shocks are in the same synchrotron regime and 
should therefore have the same color as long as p is the 
same. 

However, we do in fact observe a color change between 
the two components, of A/3 = —0.38 ±0.08 in our model. 
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There are several possible interpretations for this. 

Most likely, the forward shock is not strictly in the 
v m < v op t < v c regime at the transition time, but is in 
fact still in the process of breaking from its very blue 
rising portion (that is, v m fz v opt < v c ), and does not at- 
tain its "true" color until significantly later, when it will 
shift back to the reverse-shock color. Our data are quite 
consistent with this possibility; the physically motivated 
fit which models the reverse shock as a rising and falling 
with the appropriate change in spectral index over its 
peak is equally sound as the fit assuming no rising por- 
tion; the ultimate change in spectral index in this model 
is A@F-s,b = —0.12 ± 0.08, consistent with no change 
within 2a. If this is the case, we predict that the after- 
glow color will shift back towards the red after the end 
of our observations. 

Alternatively, one of our assumptions could be in er- 
ror. It is possible that the value of p is different between 
the forward and reverse shocks; in this case a difference 
of p rev — Pf wc i = 0.7 ± 0.1 would be required. Also, the 
amount of extinction could in fact be variable due to dust 
destruction by the GRB, though the general constancy 
of the early-time spectral index in Figure [9] before the 
transition time seems to speak against this possibility. 
Finally, if the optical or X-ray emission is strongly af- 
fected by another process such as Compton scattering, 
this could also produce spectral variability. 

Of course, a combination of two or more of these factors 
is also possible. 

4.4. Broadband Spectral Fits and Constraints on 
Extinction 

The light curve fits performed in our analysis natu- 
rally give rise to spectral fits. The filter-dependent flux 
parameters F u j give the relative fluxes in each filter for 
both components of the light curve, and can be used to 
calculate the observed SED and look for signs of host ex- 
tinction. We have already referred to the best-fit spectral 
indices (3 for the fast and slow components, all of which 
were fit assuming no host extinction. Here we will use 
different extinction models to constrain the properties of 
any host-frame dust in greater detail. 

Large amounts of extinction are implied by the X-ray 
to optical SED. In Figure [TU] we plot the SED of the 
"slow" component as computed at t = 2000 s, shortly 
before the end of our multi-filter observations and after 
the beginning of the XRT observations. 21 A model for 
Galactic extinction of E(B — V) = 0.182 mag has been 
subtracted. The predicted X-ray to optical slope, and 
the predicted slope in the optical-IR frequency window, 
is f3 = [p — l)/2 = 0.5 in the adiabatic case and also 
= 1/2 = 0.5 in the radiative case, so the prediction 
is the same. In fact, however, we measure a nearly flat 
X-ray to optical slope of Pox = 0.23 (using the i?-band 

21 Our physical model of the light curve in terms of a reverse- 
to-forward shock transition indicates that the intrinsic SED of the 
slow component may be variable due to the passage of a shallow 
minimum-energy break somewhat before this time, and so changes 
slightly after the end of our optical observations. However, to keep 
the discussion on firm observational footing without extrapolating 
our measurements beyond the range of our available data or fa- 
voring any particular model, we choose the simple empirical model 
(with no rising component) for discussion of extinction. These re- 
sults are not significantly affected by choosing the physical model 
instead. 



and 1 keV fluxes), and an IR-optical slope of f3 op t ~ 
0.95. This value of Pox is enough to unambiguously 
label this event as a "dark burst" by the criterion of 
Uakobsson et all (|2004l ) (that is, any burst with Pox less 
than the p = 2 synchrotron limit of 0.5) at this time, in 
spite of this being in fact one of the optically brightest 
bursts ever observed by Swiftl This surprising fact is 
a consequence of a combination of the unusual late-time 
X-ray brightness and the rapid early fading of the optical 
afterglow. (We could, perhaps, incorporate both this 
burst's early-time brightness and its later-time faintness 
with the moniker "grey burst," since "dark burst" seems 
inappropriate to the burst's entire evolution.) 

The most commonly invoked interpretation for these 
dark bursts, and for X-ray/o ptical flux mismatches gen- 
erally (e.g.. iSchadv et aiTl2007l ). is that the optical flux 
has been suppressed by dust extinction (or, for very high- 
redshift bursts, hydrogen absorption). A large amount of 
extinction would be necessary to apply this explanation 
here. To calculate a minimum amount of extinction, we 
assume the case that the cooling break (or for radiative 
evolution, the minimum-energy break) is just redward of 
the X-ray band, and take the minimum value of the in- 
trinsic Pox = 0.43 permitted by our uncertainty in px 
(this would imply p < 2, which is within our errors). For 
this case, the optical flux is overpredicted by more than 
a factor of 5.3 in the observed V band, requiring > 1.8 
mag of host extinction in this band, or Ay = 1 mag in 
the host-galaxy frame. This is an absolute minimum: 
requiring an intrinsic Pox = 0.5 {p = 2 or a radiative 
regime) increases this to 2.2 mag in the observed V, and 
placing the break at a more general point between the 
X-ray and optical bands increases it even further. 

Assuming a Milky Way metallicity and extinction law, 
our measurement of the host-frame hydrogen column of 
Ah, host = (1.1 ± 0.3) x 10 22 cm -2 would correspond to 
a large host-frame extinction column of Ay ~ 6 mag, 
even greater than that indicated by the X-ray/optical 
discrepancy. However, this assumption is unlikely to be 
appropriate. For a more realistic choice of the gas-to- 
dust ratio in the Small Magel lanic Cloud (SMC; about 
8.4 times the Galactic value; iGordon et al.l I2003D or a 
starburst galaxy (~ 10 times Galactic: iLisenfeld et al.1 
|2002( ). we expect a proportional reduction in Ay, pre- 
dicting Ay w 0.6-0.9 mag, roughly consistent with the 
minimum necessary suppression. 

The wide frequency coverage, photometric accuracy, 
and contemporaneous nature of our optical data gives 
us the opportunity to firmly constrain the extinction 
properties of this event. As such, we fit a power law 
plus an extinction component to our optical-IR SED, 
allowing the intrinsic unabsorbed spectral index P to 
vary as a free parameter. We try numerous extinc- 
tion pr ofiles, including the F i tzpat rick & Massa ("FM") 
model (jFitzpatrick fc Massal ll99"oT ). which has parame- 
tcrizations for a wide range of galaxy types, of which we 
fit both a Milky Way extinction pr ofile and a SMC profil e 
using the parameters measured bv I Gordon et al.l (2003). 
We also fit the Pei extinction profile |Peilll992() for the 
SMC , and the Calzetti extinction profile (jCalzetti et al.l 
2000) measured from observations of starburst galaxies. 
We use the flux parameters from our model of the slow 
component, and assume a systematic error of 2% in the 
IR filters, 3% in the optical filters, and 5% in U-band. 
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Fig. 10. — Broadband SED from optical to X-ray at t — 2000 s after the trigger. The steep X-ray decay and apparent 
shallow optical decay places the cooling break between the X-ray and optical bands. Even if we make the maximally generous 
assumption and place the cooling break at 1 keV, the optical flux is seen to be overpredicted (red line) by a factor of about 
5. If we are less generous with our assumption, and choose to interpret the late-time break in the optical afterglow seen in the 
GCN Circulars as the effect of the spectral break passing through the R band at that time, the discrepancy is even larger. The 
optical data are a good fit to a power law, and it is difficult to appeal to extinction to make up the difference. The black X-ray 
points are corrected for Galactic absorption plus a best-fit model of the host hydrogen column; grey points are corrected for 
Galactic absorption only (no host-galaxy absorption). Optical data are corrected for Galactic extinction only. 
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Fig. 11. — Plot of the minimum amount of extinction necessary in each band to resolve the discrepancy with the X-ray flux 
measurements (assuming p = 2), compared to several representative extinct i on cu rves. Ay is fixed at 1.35 mag for all models. 
The Milky Way and SMC mo dels use the formulati on of lFitzpatrick fc Massal (|1990l ) ; the Calzetti model is a model for extinction 
in starbursting galaxies from lCalzetti et al.1 (|2000f ). 
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The different extinction curves are plotted in Figure [TT1 
In every case, the best-fit model is that of no host- 
frame extinction (Ay — mag; our fits restrict Ay to 
be positive) or a very low value consistent with 0, in- 
dicating an intrinsic early-time spectral index equal to 
the observed spectral index of /3 = 0.93 ± 0.02. This 
extinction-free fit is good, with x 2 /dof = 1.91/6, and is 
plotted in Figure [P2l 

The limit on Ay depends largely on the model adopted: 
greyer dust-extinction models (and larger values of Ry, 
indicating larger grain sizes, within those models) permit 
more Ay. However, all "normal" dust-model fits strictly 
limit the observed extinction to A < 0.2 mag in the ob- 
served V band. The relatively grey Calzetti model allows 
more extinction, but even for this case, the amount of 
observed extinction can only be achieved for extremely 
large ratios of the total-to-selective extinction (Ry > 8) 
and the fit is degraded, x 2 /dof = 4.70/6 (Figure [12|) , 
though still acceptable. The \ 2 confidence contours for 
a few different fit models are plotted in Figure [T51 

This is, furthermore, for the contrived "optimistic" 
case where the X-rays are just above the break; for the 
more general case where the break is at lower energy, 
the required extinction is greatly increased. There is ev- 
idence to believe that this break is relatively near the 
optical at this time — the observations in the GCN Cir- 
culars suggest a break around t « 3 x 10 4 s, and despite 
some discrepancies all observations (including our single- 
epoch late-time PAIRITEL observations) agree that the 
optical flux by the second night has fallen well below the 
predicted value from the a < 1 decay we measure at in- 
termediate times. If we assume that this break in the 
light curve were due to a spectral break passing through 
the R band, at t = 2000 s the break would be at ^40 eV, 
predicting an optical flux 4 times (1.5 mag) higher even 
than the above "minimum" prediction would suggest, or 
3.3 mag total (Ay rs 1.8 in the host frame.) 

5. IMPLICATIONS AND ALTERNATIVE MODELS 

The standard assumptions of the fireball synchrotron 
model fare poorly when confronted with the available 
early-time data on GRB 061126. The most striking fail- 
ure is the large discrepancy between optical and X-ray 
fluxes, though the unexplained bump feature in the early 
optical light curve and the lack of agreement between 
theoretical predictions for a and (3 in X-rays are also 
unsettling. 

It is beyond the scope of this paper to fully resolve 
this quandary. Nevertheless, we present some possible 
solutions and briefly analyze their relative merits and 
failings in light of our available observations. 

5.1. Grey Dust? 

We have demonstrated that no standard extinction law 
can fully explain the discrepancy between the X-ray and 
optical fluxes at late times. Sufficiently high values of Ry 
(8-12) — corresponding to relatively grey dust laws, and 
physically to large grain sizes — do provide a respectable 
fit, but since even in the starburst galaxies studied by 
Calzetti, and in the molecular clouds in the Milky Way, 
Ry is generally in the range of only 2-5, it is worth asking 
if this is physically reasonable. 

There is ample reason to suspect that GRB host galax- 
ies, and the GRB progenitor sites within those galaxies, 



will have dust properties that do not resemble the local 
universe. GRBs in general (and GRB 061126 is no excep- 
tion) occur at high redshift and in very low-mass, high 
star-formation galaxies, and dust properties may evolve 
with time and almost certainly do correlate with host- 
galaxy type. Furthermore, if GRB progenitors are very 
massive stars, they should occur near their birth sites 
(probably within dense clouds); in the Milky Way, such 
environments are known to be correlated with anomalous 
extin ction curves and high values of Ry (|Cardelli et al.l 
1988). Finally, the GRB itself is expected to destroy 
dust within several parsecs of its progenitor, which can 
preferentially affect certain types of dust — for example, 
destroy smaller gr ains while leaving the larger grains rel- 
ative ly untouched (Waxman & Drainc 2000t iPerna et al.l 
120031 ). 

Furthermore, this is not the first GRB in which the ob- 
served optical flux has been deficient relative to the X-ray 
flux, or relative to a late-time synchrotron prediction us- 
ing the optical decay rate. In many cases a relatively 
normal extinction law [except, in nearly all cases, for the 
lack of the characteristic 2175 A bump feature seen in the 
Milky Way and the Large Magellanic Clou d (LMC)I has 
succe ssfully explained the suppression fe.g. JSchadv et all 
l2007f ) , although we caution that without broadband pho- 
tometry (including IR measurements) it can be difficult 
to accurately distinguish different models, and extinc- 
tion laws other than the standard Milky Way, LMC, and 
SMC laws are often not tested. Even so, in ma ny other 
cases (e.g. JStratta et al.ll2005UChen et al.ll2006D the low- 
rcdshift dust models clearly do not work — in these 
cases, the optical data are consistent with an unabsorbed 
power-law, but one that disagrees with an extrapolation 
from the X-ray band. Most of these authors favor a grey- 
dust interpretation of one manner or another. Unusual 
reddening laws have also been reported from observations 
of supernovae (|Wang et al1[2006l: lElias-Rosa et all [20061) 
and high-redshift quasars (|Maiolino et al.ll2004f T 

Unusual dust probably remains the simplest explana- 
tion for the flux discrepancy. Nevertheless, at least two 
factors give us pause when considering this explanation. 
First, our fits to the optical data alone simply do not re- 
quire it. There is no physically compelling reason that we 
can think of why dust should exhibit an extinction law 
that simply transforms one power law to another, requir- 
ing A\ oc (0.5 ± 0.1) \og(v/v ) + C or something mimick- 
ing it within our observational uncertainties. From the 
optical data alone, it is more natural to assume that the 
excellent power-law SED is intrinsic, and look elsewhere 
to explain the disagreement with the X-ray data. 

Second, the observed dust-to-gas ratio, while seem- 
ing to favor large amounts of extinction, actually greatly 
overpredicts the extinction if we assume a physical inter- 
pretation of the relatively grey extinction law in terms of 
dust that is mostly bound up in large grains. The grain 
size d istributions presented by IWe ingartnc r &: Draind 
(|2001h . for example, are already biased (as weighted by 
the grain volume) toward the largest grains, so for the 
same dust-to-gas ratio, skewing the grain distribution 
further toward large grains will significantly decrease the 
opacity at long wavelengths but also decrease it at short 
wavelengths. When the anticipated low dust-to-gas ratio 
of the GRB host galaxy is considered, the observed X-ray 
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Fig. 12. — Spectral energy distribution of the slowly decaying component of the afterglow at 2000 s, fit both to a model 
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Nh column is already only just enough to explain the ob- 
served extinction; suppressing it further by binding the 
available dust into very large grains becomes problem- 
atic. 

Nevertheless, without a detailed physical modeling of 
the potential dust in the GRB environment (including 
the possible effects of time-variable opacity due to dust 
destruction), we can neither firmly confirm nor rule out 
the presence of grey dust in this burst. Such an examina- 
tion is beyond the scope of our present analysis; however, 
we encourage dust modelers to make use of our observa- 
tions of the minimum extinction required (Figure Hip to 
help determine the viability of grey-dust models for this 
event. Our early-time infrared and optical photometry 
should strongly constrain any models involving dust de- 
struction. 

5.2. Synchrotron Self-Compton as the Origin of the 
X-ray Afterglow 

We now turn to interpretations in which the optical 
SED is treated as intrinsic, and instead look to processes 
primarily affecting the X-rays to explain the anomalous 
behavior. The discrepancy between the X-ray a and (3 
independently argues that the X-ray emission for this 
burst may not obey the normal assumptions about GRB 
afterglows. 

The SED in Figure [TU1 if intrinsic, has two peaks, 
reminiscent of the spectrum created by inverse-Compton 
scattering from a synchrotron source. Furthermore, 
it is interesting that the observed optical spectral in- 
dex is consistent with the X-ray spectral index: per- 
haps the X-ray afterglow is dominated by inverse syn- 
chrotron self-Compton flux boosted from the optical 
band. Compton scattering has not commonly been in- 
voked in interpreting GRB afterglow observations, but 
it has been used to explain X-ray /optical flux discrep- 
ancies similar to this one in at least two previ ously 
published cas es: GRB 000926 (lHarrison et al.ll2001h and 
GRB 030227 (jCastro-Tirado et aH\206^ ~ 

However, we consider SSC unlikely to be the solution 
to our discrepancy for several reasons. 

1. Even admitting for the effect of Compton scat- 
tering on the light curve, the optical data do not 
obey the predictions of any model without an ex- 
tinction correction. The optical and X-ray have 
very similar spectral slopes ((3 op t — 0.93 ± 0.02, 
[3 X = 1-00 ± 0.07), suggesting that the X-ray 
photons are in fact upscattered optical photons, 
and that (nongray) extinction is minimal, so this 
value of [3 must be treated as intrinsic. This re- 
stricts us to two synchrotron regimes: v c < v opt 
or v m < v op t < v c . The former case would im- 
ply p = 2 and requires a light curve that falls 
off as a = (3p - 2)/4 + (p - 2)/ (2(4 - p)) = 1.0, 
which is ruled out by the observations. (The opti- 
cal light curve decays quite slowly: a — 0.75 ±0.06 
during our observations, though it steepens later.) 
The latter case would imply p = 3, and predict 
a = 3(p — l)/4 = 1.5, which is also ruled out. 

2. Inverse-Compton scattering does not produce 
a broken power -law SED at high energies; 
ISari fc Es in (2001) have shown that in reality the 



breaks are significantly softened, in disagreement 
with the X-ray observations which indicate a sim- 
ple power-law spectral slope. 

3. Finally, it is difficult to construct the Compton- 
scattered portion of the broadband SED without 
"contaminating" the optical bands with Compton- 
scattered flux from lower energies in a way that 
would noticeably bias the blue end of our photo- 
metric SED. (However, this effect may be alleviated 
somewhat if the Compton-boosted self-absorption 
frequency v a jc is not much less than the Compton- 
boosted cooling frequency v c ,iCi since this would 
make the low-energy side of the Compton-scattered 
spectrum much steeper than what is otherwise a 
very slowly falling ((3 — —1/3) low-energy tail.) 

5.3. Physically Separate X-Ray and Optical Emission 

Regions 

Finally, we consider the possibility that the emis- 
sion regions for the optical and X-ray emission are spa- 
tially distinct. There may be two forward shocks (each 
peaking in a differ ent energy range ) resulting from a 
two-component jet (|Peng et al.ll2005l ). or the entire op- 
tical afterglow may be due to a reverse shock with 
an erratic profile extendi ng out to late times (as in 
lUhm fc Beloborodovl 120071 ) . while the X-ray light curve 
may be due to a more well-behaved radiative forward 
shock. 

The notion of finding a way to physically separate the 
optical and X-ray emission has recently found support 
in other contexts, in particular to try to resolve the ap- 
parent lack of correlation betw een Swift X-ray and op- 
tical breaks (jOates et al.ll2007f) and to explain Swift X- 
ray breaks generally (|Panaitescull2006[ ). Unfortunately, 
the relatively short period of overlap between the optical 
and X-ray data for GRB 061126 prevents us from being 
able to constrain this possibility in any detail for this 
event. However, we will note that any two-origin expla- 
nation shares the difficulty just discussed in the context 
of whether inverse Compton may explain the X-ray flux 
excess — that is, how to avoid contaminating the blue 
end of the observed optical spectrum with emission from 
the low-energy tail of a synchrotron-like spectrum peak- 
ing at X-ray wavelengths. 

6. CONCLUSIONS 

We have presented multi- wavelength, early-time obser- 
vations of the recent, bright GRB 061126. In our favored 
model, the early optical data appear to be well explained 
by a reverse shock. The reverse-shock emission fades 
rapidly and is dominated by emission from the forward 
shock from ^500 s onward. A small color change is ob- 
served at this transition, likely due to the peak of the 
forward-shock synchrotron spectrum passing through the 
optical band at roughly the same time. 

However, only with significant contrivance can this 
model also explain the X-ray observations from the Swift 
XRT, which are much brighter than an extrapolation of 
the optical flux would predict. The optical flux may be 
suppressed by host-galaxy extinction, but this would re- 
quire a remarkably grey dust law, capable of generat- 
ing 1.5-3 mag of optical extinction without causing the 
observed U-K s SED to deviate from a power law. Al- 
ternatively, the X-ray and optical emission may be due 
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to two physically or spatially distinct emission processes, 
but our observations strongly restrict the positions and 
sharpness of any spectral breaks and cast doubt on the 
viability of this model. 

In either case, we advise caution in interpreting future 
early afterglows in the absence of data sets that are not 
well sampled both spectrally and temporally. Many pre- 
vious studies, by necessity, have been restricted in cither 
their temporal properties (sampling) , frequency coverage 
(available filters and wavebands), or both. Under the 
standard fireball model these assumptions seemed valid, 
as only a single parameter set would fit them. Here we 
show that when a more complete picture is available, no 
parameter set seems to fit the data, unless we invoke 
large quantities of grey dust or separate origins for the 
X-ray and optical afterglows. 

The very bright X-ray afterglow and late-time optical 
faintness are enough to qualify this event under some 
definitions as a dark gamma-ray burst in spite of the ex- 
traordinary bright early afterglow. Had the event been 
intrinsically fainter or at higher redshift, or had the opti- 
cal follow-up observations been delayed significantly, the 
optical afterglow may have been missed entirely. There- 
fore, it is possible that events like GRB 061126 may rep- 
resent a significant fraction of dark bursts. 

The most common interpretation of dark bursts is that 
they are due to absorption of the optical flux by host- 
galaxy dust or by neutral hydrogen at very high redshift. 
Our study presents tentative evidence for a third pos- 
sibility, which is that the optical faintness may be in- 
trinsic to the GRB itself, due to enhancement of the X- 
ray flux or intrinsic suppression of optical flux, or both. 
If the optical faintness is due to dust, our results sug- 
gest that near-IR and optical observations alone may 
not be as constraining as once hoped, since any extinc- 
tion must not cause significant deviations from a power 
law, even across a decade in frequency from the near- 
IR to the UV. However, the combination of near-IR, 
optical, and X-ray observations remains a potent tool 
for understanding GRBs, and we anticipate that addi- 
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TABLE 6 
GCN Observations of GRB 061126 



GCN 




Filter 


Magnitude 13 


Flux b 


Re 




(s) 




(MJy) 




5869 


36.3 


R 


12.93 ±0.20 


20893 ± 3515 


(1) 


5857 


258.3 


R 


15.97 ± 0.05 


1270.6 ± 57.2 


(2) 


5868 


626.4 


1 


17.10 ± 0.20 


351.6 ± 59.1 


(3) 


5868 


806.1 


1 


17.00 ± 0.20 


385.5 ± 64.9 


0) 


5859 


2820.1 


r> 

H 


1 o on _l_ n on 
18.80 ± 0.20 


CiQ O _1_ 1 K Q 

93.8 ± 15.8 


/ i \ 

(v 


5859 


5880.4 


R 


19.30 ±0.50 


59.2 ±21.8 


("1) 


5859 


8939.8 


R 


19.40 ± 0.30 


54.0 ±13.0 


("1) 


5866 


34367.3 


R 


21.50 ±0.25 


7.8 ± 1.6 


(5) 


5876 


152064 


R 


22.85 ±0.06 


2.2 ±0.1 


(6) 


5875 


156381 


R 


23.69 ±0.17 


1.0 ±0.2 


(7) 


5902 


39225.6 


R 


21.16 ±0.04 


10.7 ±0.4 


(«) 


5902 


51235.2 


R 


21.65 ±0.08 


6.8 ±0.5 


(8) 


5903 


35424.0 


R 


20.98 ±0.10 


12.6 ±1.1 


(9) 


5903 


36288.0 


R 


21.04 ±0.09 


11.9 ±0.9 


( n ) 


5903 


42336.0 


R 


21.34 ±0.10 


9.0 ±0.8 


(9) 


5903 


48384.0 


R 


21.49 ±0.10 


7.9 ±0.7 


(9) 



Note. - Ref eren ces. - (1) I Williams fc Milnd 20061: (2) 
Smith et all [20061 ; (3) iToriil [20061; (4) lUpdike etldl p006l; (5) 
Kannfc Malesanil 120061 ; (6 ) iRol et at] 120061 ; (7) iKann 200j; (8) 



Pozanenko et al.ll2006l ; (9) lMisral l2006 



Exposure mid-time, measured from the Swift trigger (UT 
08:47:56. 4). b Observed value; not corrected for Galactic extinc- 
tion. 



TABLE 7 

Results of Band et al. (1993) model fits to the BAT+RHESSI spectrum of 

GRB 061126 



Region 


a 


E oha , 

" peak 

[keV] 


100 keV Norm, 
[ph cm -2 s _1 ] 


Xl (DOF) 


Full (t = -6-35s) 
Pulse A (t = 3-14s) 
Pulse B (t = 19-25s) 


-1.06 ± 0.07 

-0.94 ±0.06 

-0 9+ - 2 
u - 9 -o.i 


< -2.3 620tig° 

< -2.5 790+130 

350t»° 


7.8l°'5 x 10~ 3 
1.8 ±0.1 x 10" 2 
l.l +0 i x 10- 2 


0.733 (106) 
1.042 (105) 
1.179 (100) 



Note. — The quoted uncertainties correspond to the 90% confidence region. The 
data in each time region are acceptably fit by an exponential times a power-law 
model. The high-energy power-law component (with photon index 0) is not required 
in the fits but can be constrained for regions "Full" and "Pulse A." Using a m — 1 and 
the declining i? pca k and normalization values between pulses A and B, we estimate 
an approximately energy-independent GRB spectral flux of 0.1+ 'g mjy below 1 keV 
at t = 60 s (Figure 0. 



TABLE 8 
Optical Light-Curve Fits 







C*S,6 a S,a 




X 2 /dof 


2.09 ±0.29 


1.28 ±0.01 


0.80 ±0.05 


-0.32±0.03 a 


212.0/185 


1.76 ±0.22 


1.07 ±0.02 


0.58 ±0.12 


b 


261.8/193 


1.96 ±0.29 


1.31 ±0.02 


0.75 ± 0.06 


-0.38 ± 0.08 


231.3/192 


1.70 ±0.09 


1.23 ±0.02 


-0.50 b 0.75 b 


-0.12±0.09 c 


227.6/197 



Note. — Summary of key parameters and \ 2 from the various models 
fit to the data. The model with no color change is strongly ruled out. 
The nature of the color change depends on the assumed model. 
a Not a formal fit parameter in this model — the flux amplitude parame- 
ters in each filter are allowed to assume their arbitrary best-fit values. In 
other models the change in these parameters is constrained to be due to 
variation in the spectral index j3. b Fixed parameter. Change between the 
spectral index of the fast component and the index fit to the slow com- 
ponent after its peak (/3s b)- The slow component undergoes a chromatic 
break from @ s .b = 0.284 'to p s ,b = 1-11. 



